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Abstract 

The spin asymmetry parameter A T characterizing the angular distribution of the 
total hadron momentum in the decay of a polarized tau can be calculated rigorously 
using perturbative QCD and the operator product expansion. Perturbative QCD cor- 
rections to the free quark result A T = 1/3 can be expressed as a power series in a s (M T ) 
and nonperturbative QCD corrections can be expanded systematically in powers of 
1/M?. The QCD prediction is A T = 0.41 ± 0.02. In the decay of a high energy tau 
into hadrons, the value of the hadronic energy distribution dR T /dz evaluated at the 
maximum hadronic energy fraction z = 1 can also be calculated rigorously from QCD. 



The spin-dependence of processes involving elementary particles contains a wealth 
of information about their fundamental interactions. Unfortunately this information is not 
easily accessible to experiment. It requires either the use of polarized beams and targets, or 
the measurement of the polarization of final state particles. The tau lepton is one of the few 
elementary particles whose spin can be effectively analyzed by its decays. It is well known 
that several of the exclusive decay modes of the tau can be used to analyze its spin [||, |2| . In 
this Letter, I point out that inclusive decays into hadrons can also be used for this purpose. 
The polarization of a sample of taus that decay into hadrons can be determined from the 
angular distribution of the total momentum of the hadrons. The asymmetry parameter that 
characterizes the angular distribution can be computed systematically using perturbative 
QCD and the operator product expansion. In a sample of high energy taus that decay into 
hadrons, the endpoint value of the hadronic energy distribution can also be computed rigor- 
ously from QCD, and can be used to measure the average helicity of the taus. Measurements 
of these spin-dependent observables could provide dramatic confirmation of the applicability 
of perturbative QCD to the inclusive hadronic decays of the tau. 

It is convenient to normalize the inclusive decay rate of the tau lepton into a neutrino 
plus hadrons to the electronic decay rate by defining the ratio 



For a sample of taus with polarization P, the angle 6 between the total momentum of the 
hadrons in the tau rest frame and the spin quantization axis has a distribution proportional 
to 1 + A T P cos 9, where A T is an asymmetry parameter. This angular distribution can be 
used to separate R T into "forward" and "backward" components Rf and Rb'- 
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The asymmetry parameter A T is then 



Rf + Rb 



(3) 



A naive estimate of the asymmetry parameter can be obtained by considering the 
decay of the tau into hadrons at the quark level, where it proceeds through the processes 
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t~ — > v T du and r~ — > v T su. The momentum of the du and su pairs can be identified with 
the total momentum of the hadrons. Ignoring the QCD interactions that bind the quarks 
into color singlet hadrons, the angular distribution of the total hadron momentum is 



^(|Kd| 2 + |Ks| 2 ) 1 + -Pcos9) . (4) 
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The squares of the Kobayashi-Maskawa matrix elements add up to 1 to high accuracy, so they 
will be omitted below. The naive estimates for the ratio ([!]) and the asymmetry parameter 
(D are R T = 3 and A T = 1/3. 

In the case of R T , the QCD corrections to the naive result can be computed systemati- 
cally using perturbative QCD and the operator product expansion || £|, |[. The perturbative 
corrections can be expanded as a power series in a s (M T ) |§ and the nonperturbative cor- 
rections can be organized systematically into an expansion in powers of 1/Mj?. A thorough 
analysis of the QCD and electroweak corrections to the ratio R T has recently been carried 
out j7|. The methods that were used to calculate the ratio R T can also be used for a rigorous 
calculation of the asymmetry parameter A T . The starting point is an expression for the 
angular distribution of the total hadron momentum as an integral over the invariant mass s 
of the hadrons: 



dR T 
dcosO 
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+ Imn (0) (s + ie)(l + P cost/ , 



(5) 



where n^ J ^(s), J = 0, 1 are the transverse and longitudinal correlators for the quark current 
that couples to the virtual W boson. The notation is the same as in Ref. [0]. The correlators 
n^ J ^(s) are analytic functions of s except along the positive real s-axis. This allows the 
integral in (|5|) to be expressed as a contour integral in the complex s-plane. The contour 
can be deformed so that it runs counterclockwise around the circle \s\ = M^. The resulting 
expressions for the forward and backward components of R T defined in (0) are 
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The contour integral expressions (@) and ([?D reveal that the polarization asymmetry A T , like 
the ratio R T , is completely determined by correlation functions at the distance scale 1/M T . 
This implies that the nonperturbative long distance effects of QCD can be expressed in terms 
of matrix elements of local operators. These matrix elements appear when the operator 
product expansion is used to expand the correlators n( J )(s) in @ and ([?]) in powers of 1/s. 
Evaluating the contour integrals, the QCD corrections to the naive predictions Rp = 2 and 
Rb = 1 are obtained as systematic expansions in powers of There is also an important 

electroweak correction consisting of a multiplicative short distance factor Sew — 1-019 ||. 
The resulting expressions for the forward and backward components of R T have the form 

R F = 2S EW (l + SP + # + 4 4) + 4 6) + •••) , (8) 

r b = s EW (i + 4 0) + 4 2) + 4 4) + 4 6) + •••) , (9) 

where the fractional corrections 8^ and 8^ are proportional to 1/M™ with coefficients that 
depend logarithmically on M T . For R T = Rp + Rb, the fractional corrections to the free 

quark value 3 Sew are (25^ + 8^)/3. 

The fractional corrections 8p and can be calculated straightforwardly using the 
operator product expansions for the correlators IB J )(s) that are collected in Ref. 0. The 
dimension-0 corrections, which represent the purely perturbative effects from the interactions 
of massless quarks and gluons, are 

4 0) = ^ + 5.765 (^) 2 + 34.48 (^Y + (d 4 + 165.1) f^V , (10) 

4 0) = ^ + 4 - 077 (v) 2 + 10 - 13 (^) 3 + ( rf 4-96.1)(^) 4 , (11) 

where a s = a s (M r ) is the running coupling constant of QCD in the MS scheme evaluated at 
the scale M T . The coefficient d 4 in the term in (|i~0l) and (|TTD is the fourth coefficient in the 
perturbative expansion of — 2n 2 s(d/ ds)!!^^) in powers of a s /n and has not been calculated. 
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The previous coefficients are 1, 1, 1.64, and 6.37 ||. We assign a very conservative error 
to this unknown coefficient: c/4 = ± 100. The corresponding coefficient in the fractional 
correction to R T is then 78.0 ± 100. The dimension-2 corrections are perturbative corrections 
due to the running quark masses. The only correction that is numerically significant comes 
from the strange quark mass m s = m s (M T ): 



5y = -6 sin 1+ — 2 , (13) 



3 W M r 2 

where #c is the Cabbibo mixing angle: sin 2 9 C = 0.049. For the running strange quark mass 



in the MS scheme evaluated at the scale M T , we take the value m s (M T ) = (0.17±0.02) GeV. 
The first nonperturbative corrections appear at dimension 4 in the form of scale invariant 
matrix elements called the gluon condensate and quark condensates: 

4" = 2^ (l - "g) + 48^^ - gUfeJl £ , (14) 



18W M r 4 7 " a s M T 4 

4 4) = -4^ (1 - 11^ <( " s/ ;j GG> . as) 
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Note that the contribution of the gluon condensate <(a s /ir)GG> cancels to order a s in the 

fractional correction (28 p +Sg )/3 to the ratio i? T . This results in a suppression of the gluon 
condensate contribution to R T by two orders of magnitude. We take the value of the gluon 
condensate to be <(a s /n)GG>= (2 ± 1) x 10~ 2 GeV 4 [10|. The matrix element <mipip> in 
( p!4|) is a weighted average of the quark condensates: 

Tl <m u uu> + cos 2 9c <mddd> + sin 2 6c <m s ss> , s 

<mipip> = . (16) 

Its value is <rm]j?p>= (—8 ± 1) x 10 -5 GeV 4 . The inverse power of a s (M T ) multiplying the 
m s (M T ) A term in ( |l4|) was first understood by Broadhurst and Generalis [|TTJ. At dimension 
6, there are too many unknown matrix elements for a completely systematic treatment. 
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Within the vacuum saturation approximation, these corrections are 



(6) 256tt 3 pa s 2 f . 

8f = m*—, (17) 

(6) 20487T 3 pa s <^> 2 
S B ~ ' (18) 

with p — 1. It has been found empirically that this approximation underestimates the size of 
the dimension-6 correction, and it is better to treat pa s <ipip> 2 as an effective scale-invariant 
operator of dimension 6, independent of <ipip>. The best estimate for this parameter is 
pa s <ipip> 2 = (4 ± 2) x 10~ 4 GeV 6 . The dimension-8 and higher corrections are assumed to 
be completely negligible. 

Inserting the fractional corrections given above into (||) and we obtain predictions 
for the ratio R T = Rp + Rb and the asymmetry parameter A T defined in (Q) as a function 
of a s (M T ) and the five parameters d^, m s (M T ), <(a s /7i)GG>, <mipip> ) and pa s <iptfj> 2 . 
Alternatively, given a value for R T , we can predict both a s (M T ) and A T . The predictions are 
shown in Table 1. The uncertainty in a s (M T ) is dominated by the assumed error of ±100 
in the coefficient <i 4 . For R T = 3.60, the uncertainty in a s (M T ) is 3.6%. After <i 4 , the next 
largest errors are 1.1% from pa s <i\) r i\)> 2 and 0.4% from m s (M T ). The uncertainty in A T is 
dominated by the gluon condensate, and is 5.4% for R T = 3.60. The next largest errors are 
1.8% from pa s <ipip> 2 and 0.7% from d^. 

There are two independent ways of measuring the ratio R T experimentally. Using the 
universality of electron and muon couplings, it can be expressed in terms of the electronic 
branching fraction B e of the tau: R T = (1/B e ) — 1.973. Alternatively, using the universality 
of electron and tau couplings as well, it can be expressed in terms of the masses and lifetimes 
of the /i and r: R T = (t^/t t )(M^/ M T ) 5 — 1.973. The present world average for the electronic 
branching fraction is B e = (17.78±0.15)% |T2|, and it gives the ratio R T = 3.651 ±0.047. The 
present world average for the tau lifetime is r T = (2.96 ±0.03) x 10~ 13 s [12]. Combined with 



the recent precise measurement of the tau mass jCJ, it gives the ratio R T = 3.545 ± 0.056. 



Forming the weighted average of the two independent determinations of R T , we get R T = 
3.607 ± 0.036. From Table 1, this determines the running coupling constant at the scale 



M T to be a s (M T ) = 0.319 ± 0.017. We have added in quadrature the error from Table 1 
and the error due to the experimental uncertainty in R T . Evolving the running coupling 
constant up to the Z° mass, it becomes a s {Mz) — 0.1176 ± 0.0021. The QCD prediction for 
the asymmetry parameter is A T = 0.413 ± 0.022. A measurement of A T consistent with this 
prediction would provide dramatic confirmation of the accuracy of the QCD calculation of 
the ratio R r . 

The angular distribution dR T /dcos9 that defines the asymmetry parameter A T is 
most easily measured in low energy experiments. In high energy experiments, such as the 
decay of the Z° into taus, the quantity that is most easily measured is the distribution 
dR T /dz of the hadron energy fraction z = Eh/E t , where Eh is the total energy of the 
hadrons in the rest frame of the Z° and E T is the energy of the decaying tau. If the spin 
quantization axis is chosen to lie along the direction of the tau momentum, then the variables 
z and cos 6 are related by 

(2z - l)Ml - s 

for a hadronic state with invariant mass s [0. Changing variables in (|5|) from cos 9 to z, we 
obtain an expression for the hadronic energy distribution Dp(z) = dR T /dz for a tau with 
polarization P in the helicity basis. At the endpoint z — 1, it reduces to 



D P (z = i) = iar jf A _ _L) (i mn « (s + fc) (i + p + 2^(1 - P) 



+ Imn (0) (s + k) (1 + P)^j 



(20) 



We define the energy asymmetry function A(z) by 

D P (z) = D (z) (1 + A(z)P) , (21) 

where D (z) is the energy distribution for an unpolarized tau. 

By the same arguments that were used for R T and A T , the quantities -D (l) and 
A(l) can be calculated rigorously using QCD. Perturbative corrections can be expanded as 
a power series in a s (M T ), and nonperturbative corrections can be expanded systematically 
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in powers of 1/M 2 . The free quark prediction are -D (l) = 5 and A(l) = 1/5. We express 
the QCD corrections in terms of fractional corrections 5^ and 8^ which scale like 1/M": 

D P=+1 (z = 1) = 6 S EW (l + 4 0) + 4 2) + 4 4) + 4 6) + ...) , (22) 

^-1(2 = 1) = 4 S w (l + 4 0) + 4 2) + 4 4) + 4 6) + •••) • (23) 

The fractional corrections due to the perturbative interactions of massless quarks and gluons 
are 



V + 5 - 015 (v) +2450 (v) +^ + 68.7)(^) , (24) 
4 0) = ^ + 3.515 + 4.54 + (d 4 - 123.9) . (25) 



The fractional corrections from the running strange quark mass are 

13 a. 



4, = _ 6sin ^( 1 + ^)- (26) 
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At the order in a s (M T ) to which we have calculated the coefficients of the matrix elements, 
the fractional corrections 8$ and 8^£\ n — 4, 6, are related in a simple way to the fractional 
corrections to Rp and Rb- The dimension-4 corrections are 8$ = 5^/3 and 5^ 4 ' ) = 8$ /2 and 

the dimension-6 corrections are 5^ = and 5^ = 8$/ 4. The QCD predictions for D (l) 
and A(l) as a function of R T are presented in Table 1. The error on D (l) is remarkably 
small, only 0.4% for R T = 3.60. The dominant errors are 0.3% from <(a s /7r)GG>, 0.3% from 
g?4, and 0.1% from pa s <tj)ip> 2 . The error in A(l) is dominated by the gluon condensate and 

is 4.5% for R T = 3.60. The next largest errors are 1.1% from d A and 0.9% from pa s <tjj^> 2 . 
Taking R T = 3.607 ± 0.036 and combining the error from Table 1 in quadrature with the 
error due to R T , the QCD predictions are D (l) = 5.853 ± 0.051 and A(l) = 0.246 ± 0.011. 

The QCD predictions for A T and for the high energy parameters D (l) and A(l) are 
sufficiently precise that measurements of these parameters could be used to determine the 



polarization of a sample of taus. The polarization of the taus produced in the decay of the 
Z° has been measured using the energy distributions of the exclusive decay modes e~u e u T , 
^u^t, 7r~z/ r , p~z/ r , and a{v T [EH]. The measurement of the endpoint value Dq(1)(1 + 
A(1)P) of the inclusive hadronic energy distribution might be competitive as a method 
for determining the tau polarization P, since it is less sensitive to errors due to particle 
identification. 

Precise measurements of the asymmetry parameters could also be used to determine 
the nonperturbative matrix elements that arise in the operator product expansion. The 
errors in A T and ^4(1) are dominated by the gluon condensate <(a s /7r)GG>, with the other 
errors being smaller by a factor of 3. Thus precision measurements of A T and A(l) could 
improve the determination of the gluon condensate by about a factor of 3. 

In this Letter, we have introduced several new observables involving hadronic decays 
of the tau lepton that can be calculated rigorously from QCD. The only ones that were 
known previously were the ratio R T and the moments of the invariant mass distribution 
dR T /ds j|, nj. The endpoint value of the hadronic energy distribution dR T /dz in the decay 
of a high energy tau can also be calculated rigorously, as can the moments of the energy 



distribution |L6|. For each of these quantities, there is also a polarization asymmetry that 
can be calculated. The QCD predictions for the quantities A T , D (l), and -4(1) differ from 
the free quark predictions by about 20%, which is much larger than the errors in the QCD 
predictions. Thus hadronic decays of the tau lepton provide a remarkable laboratory in 
which QCD can be tested at low energies with unprecedented precision. 

This work was supported in part by the U.S. Department of Energy, Division of High 
Energy Physics, under Grant DE-FG02-91-ER40684. 
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Table Captions 

l. QCD predictions for a s (M T ), A T , D (l), and A(l) as a function of the ratio R T . The 
errors due to variations of d 4 , m s (M T ), <(a s /n)GG>, <mijji/j> , and pa s <ipip> 2 have 
been added in quadrature. 



R T 


a s {M T ) 


A T 


A>(1) 


A(l) 


3.50 


0.287 ±0.009 


0.407 ±0.023 


5.735 ±0.022 


0.241 ±0.011 


3.52 


0.294 ±0.009 


0.408 ± 0.023 


5.761 ±0.022 


0.242 ±0.011 


3.54 


0.301 ±0.010 


0.409 ± 0.023 


5.786 ±0.023 


0.243 ±0.011 


3.56 


0.308 ±0.010 


0.411 ±0.023 


5.812 ±0.023 


0.245 ±0.011 


3.58 


0.315 ±0.011 


0.412 ±0.023 


5.837 ±0.024 


0.246 ±0.011 


3.60 


0.321 ±0.012 


0.413 ±0.022 


5.862 ±0.025 


0.247 ±0.011 


3.62 


0.328 ±0.012 


0.414 ±0.022 


5.886 ±0.026 


0.248 ±0.011 


3.64 


0.334 ±0.013 


0.415 ±0.022 


5.911 ±0.027 


0.250 ±0.011 


3.66 


0.340 ±0.013 


0.417 ±0.022 


5.935 ±0.028 


0.251 ±0.011 


3.68 


0.346 ±0.014 


0.418 ±0.022 


5.959 ±0.029 


0.252 ±0.011 


3.70 


0.352 ±0.015 


0.419 ±0.022 


5.983 ±0.031 


0.254 ±0.011 



Table 1 
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